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ABSTRACT 

Dynamical systems in equilibrium have a stationary entropy; we suggest that ellipti- 
cal galaxies, as stellar systems in a stage of quasi-equilibrium, may have in principle 
a unique specific entropy. This uniqueness, a priori unknown, should be reflected in 
correlations between the fundamental parameters describing the mass (light) distribu- 
tion in galaxies. Following recent photometrical work on elliptical galaxies (Caon et 
al. 1993; Graham & CoUess 1997; Prugniel & Simien 1997), we use the Sersic law to 
describe the light profile and an analytical approximation to its three dimensional de- 
projection. The specific entropy is then calculated supposing that the galaxy behaves 
as a spherical, isotropic, one-component system in hydrostatic equilibrium, obeying 
the ideal gas state equations. We predict a relation between the 3 parameters of the 
Sersic law linked to the specific entropy, defining a surface in the parameter space, 
an 'Entropic Plane', by analogy with the well-known Fundamental Plane. We have 
analysed elliptical galaxies in two rich clusters of galaxies (Coma and ABCG 85) and 
a group of galaxies (associated with NGC 4839, near Coma). We show that, for a 
given cluster, the galaxies follow closely a relation predicted by the constant specific 
entropy hypothesis with a typical dispersion (one standard deviation) of 9.5% around 
the mean value of the specific entropy. Moreover, assuming that the specific entropy is 
also the same for galaxies of different clusters, we are able to derive relative distances 
between Coma, ABGC 85, and the group of NGC 4839. If the errors are only due to 
the determination of the specific entropy (about 10%), then the error in the relative 
distance determination should be less than 20% for rich clusters. We suggest that the 
unique specific entropy may provide a physical explanation for the distance indicators 
based on the Sersic profile put forward by Young & Currie (1994, 1995) and recently 
discussed by Binggeli & Jerjen (1998). 

Key v^rords: galaxies: clusters: individual: Coma, ABCG 85 - galaxies: distances and 
rcdshifts - galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters - 
distance scale. 



1 INTRODUCTION 

Elliptical galaxies are believed to be self-gravitating systems 
in a quasi-equilibrium state. They are far from being the 
simple, dull objects they were thought to be until the early 
70's: high spatial and spectroscopic observations have re- 
vealed that elliptical galaxies show a great variety of fine 
structures (e.g. Kormendy 1984, Bender & MoUenhoff 1987) . 
However, underlying these fine structures, elliptical galaxies 
have a striking regularity concerning their radial light distri- 
bution. Indeed, until the mid 80's, it seemed that the bright- 



ness profile of elliptical galaxies could be well described by 
the de Vaucouleurs (1948) R^^"^ law (e.g. Kormendy 1977). 
Nevertheless, systematic departures from the R^^* law were 
observed (Michard 1985; Schombert 1986, 1988) and, more 
recently, thanks to higher quality observations, it was shown 
(Caon et al. 1993; Graham & CoUess 1997) that elliptical 
galaxies in a very wide range of sizes can be well modelled 
by the Sersic (1968) empirical law (expressed in magnitude): 

M = MO + fcn(^)'^/"^ (1) 
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a non-homologous generalization of the R}'^'^ law. is the 
central brightness, ReS is the effe ctive radius and fe„ a func- 
tion depending on n (see below, §3.1). 

Besides the light profile, another very important reg- 
ularity observed in elliptical galaxies is the so-called Fun- 
damental Plane (Dressier et al. 1987b; Djorgovski & Davis 
1987). Elliptical galaxies populate a two dimensional sur- 
face in the three dimensional space defined by the effective 
radius, RcS, the mean surface brightness inside this radius, 
(/i)off, and the central line of sight velocity dispersion, cro. 
The understanding of such a relation should provide impor- 
tant clues to the formation and evolution of these galaxies. 

What is the origin of the relations among the parame- 
ters describing the light (and mass) distributions of elliptical 
galaxies? If these objects are, as we believe, in equilibrium, 
then the virial theorem should apply and we can derive for 
a given mass distribution a relation among the gravitational 
radius, total mass, and mean velocity dispersion. This rela- 
tion can be 'translated' in the observed Fundamental Plane 
with the help of some hypotheses related to the mass to 
light ratio, light distribution, and hydrostatic equilibrium 
(Bender et al. 1992, Ciotti & Lanzoni 1997). The validity of 
these hypotheses is, however, still a matter of debate (cf. for 
instance Prugniel & Simien 1997). 

Other than its theoretical importance, the scaling rela- 
tions of the fundamental parameters of galaxies may have a 
practical application as a distance indicator. Thus, like the 
Fundamental Plane (Djorgovski & Davis 1987), the Dn-cr re- 
lation (Dressier et al. 1987a), Tully & Fisher (1977) method 
for spiral galaxies, etc., the parameters of the Sersic law 
can also be used to derive distances (Young & Currie 1994, 
1995). 

In this paper we will look for fundamental relations 
among the parameters describing the light distribution of 
elliptical galaxies. The starting point is also the fact that 
these galaxies are in equilibrium, but we will use the max- 
imum entropy postulate of thermodynamics instead of the 
virial theorem. We will derive an expression relating the free 
parameters of the Sersic profile to the specific entropy. We 
will also give a number of useful and simple expressions to 
calculate the astrophysical quantities related to the Sersic 
profile and its three dimensional deprojection. We will then 
show that the specific entropy relation holds for at least two 
clusters of galaxies and we suggest that it may, in principle, 
be used as a relative distance indicator like the shape-profile 
distance indicators introduced by Young & Currie (1994, 
1995). 



2 THE ENTROPY 
2.1 Thermodynamics 

The regularity observed in the light distribution of elliptical 
galaxies may be conveniently understood as an evolutionary 
convergence. Independently of the initial conditions, these 
galaxies evolve towards a single class of objects. Since the 
stellar two-body relaxation time scale is far greater than 
the age of the Universe, the mechanism behind this evolu- 
tion is probably the so-called violent relaxation (Lynden-Bell 
1967). The rapidly changing global gravitational potential 
(due to an initial collapse or merging) produces a mixing in 



phase space within the same time scale as the collapse or 
merging itself. 

The second postulate of thermodynamics says that, 
when a dynamical system is in equilibrium, then this system 
is in a state of maximum entropy. If, for a given system, there 
exists only one state of thermodynamic equilibrium that can 
be reached, independently of the initial conditions, then it 
is associated with a stable state of single, global maximum 
entropy. Gravitational systems, however, present some dif- 
ficulties: given an arbitrary isolated, self-gravitating system 
of point particles, the entropy can grow indefinitely by col- 
lapsing part of the mass in the centre and expanding the 
outer parts to infinity, while still conserving the total mass 
and energy. This is the gravo-thermal catastrophe which can 
actually occur because a gravitational system has negative 
specific heat (Lynden-Bell & Wood 1968; Saslaw 1985). 

Thus, if we try to maximise the entropy (requiring that 
the mass and energy be constant) in order to get the equi- 
librium (and hopefully stable) configuration, we obtain an 
isothermal sphere (Lynden-Bell & Wood 1968), which does 
not match the observed elliptical galaxies. A lot of work has 
been done these last three decades to obtain the structure 
of self-gravitating systems from the maximum entropy prin- 
ciple (e.g. Antonov 1962; Lynden-Bell 1967; Shu 1978; Bin- 
ney 1982; Tremaine et al. 1986; Dejonghe 1987; Richstone 
& Tremaine 1988; Merritt et al. 1989; Spergel & Hernquist 
1992; Hjorth & Madsen 1995). A local maximum entropy can 
be found if some ad hoc constraints are imposed (e.g. im- 
posing a power law density profile. White & Narayan 1987). 

Here, we adopt a different approach: we just assume the 
existence of a local maximum entropy associated with the el- 
liptical galaxies, without trying do obtain this entropy from 
first principles. Except for interacting galaxies, it is reason- 
able to suppose that elliptical galaxies are, at least, objects 
in thermodynamic quasi-equilibrium (the best evidence be- 
ing the regularity in the luminosity distribution) . Hence, we 
will proceed computing the entropy of a self-gravitating sys- 
tem using classical thermodynamics. Let us assume that the 
stars in a galaxy behave like an isolated, self-gravitating gas 
of equal mass particles. Consequently, we start by taking the 
differential of the entropy: 



dS - 



LdE+-dV-^dN, 



(2) 



where S, E, V, N, T, P, fiq are, respectively, the entropy, en- 
ergy, volume, number of particles, temperature, pressure, 
and chemical potential. Formula (b|) is nothing more than 
the fundamental equation of thermodynamics (also know 
as the Euler relation) in the entropy representation (e.g. 
Callen 1960), which applies to any thermodynamic system 
in equilibrium. For an isolated, self-gravitating system, the 
extensive and intensive parameters in Eq. (Q) are usually 
interpreted using their classical meaning. For instance, the 
chemical potential fiq (that appears already in the work of 
Lynden-Bell 1967 in the context of stellar dynamics) gives 
the variation rate of the energy with the number of particles, 
di5/dA'^, for constant entropy and volume (see also Saslaw 
1985) or, from the statistical physics point of view, fiq is re- 
lated to the Fermi energy of the system (Stiavelli 1998). The 
temperature is defined in terms of the velocity dispersion a, 
as kT = ma^ (m is the mean particle mass of the system). 
We now need the equations of state in order to describe 
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the properties of tlie 'gas' of stars. Tliere is no known exact 
equation of state for a gas of self-gravitating particles, so 
we will use the ideal gas equations which are usually a good 
approximation for hot, structureless, rarefied gases: 

3 



E 



-N kT : 



feTln 



PV = N kT : 





3/2\ 


2nm kT 





(3) 



where k and h are the Boltzmann and Planck constants. 
Although the above equations of state are a good approx- 
imation for non-mteracUng, mono-atomic, hot gas (in fact, 
these equations violate the third law of thermodynamics or 
Nernst postulate, i.e. that the S —> as T 0), their 
application for interacting particles is doubtful. Modifica- 
tions of the Eqs. have been proposed; Bonnor (1956, and 
references therein) discusses a modification of Boyle's Law 
that takes self-gravity into account. However, that modifi- 
cation was derived for a bounded, isothermal sphere and 
it adds considerable complication to the problem without 
really putting away the simplifying hypotheses. Therefore, 
we take a pragmatic approach and we will use the simplest 
physically acceptable equations of state, Eqs. (^. The va- 
lidity or not of our approach should be verified a posteriori. 

Thus, we can now write the entropy as a function of the 
mass density p = mN/V and pressure P. Using Eq. ^ and 
Eqs. (^, we obtain, after straightforward algebra: 

We now define the specific entropy s = S/M, where M is 
the total mass of the system. Therefore we write: 



-5/2 p3/2\ 



p dV + constant , 



(5) 



where the 'constant' above does not depend on the density 
distribution and we will drop it hereafter. Notice that the 
above equation is the same used by White & Narayan (1987) . 
Equation (^) can be written in a more compact form using 
the definition of a mean value, (x) = J xf{x)dx/ J f{x)dx, 
and interpreting the density p as a distribution function: 



(ln(, 



-5/2 p3/2 



)) 



(6) 



In Appendix 1 we show that we can obtain the same form for 
the specific entropy as in Eq. (^ using the Boltzmann- Gibbs 
definition of the entropy. 



2.2 Entropy and light profile 

The parameters needed for the evaluation of the specific 
entropy of early type galaxies are not accessible by direct 
observation. The calculation of the specific entropy must 
thus be obtained using the available light profile and some 
additional hypotheses. 

It is clear that, to calculate the mass distribution needed 
to compute the specific entropy of the galaxy from the ob- 
served light distribution, we will need to introduce a mass- 
to-luminosity ratio, M/L. In general, the radial variation 
of the M/L ratio in elliptical galaxies is poorly known, ba- 
sically because its determination depends on knowledge of 
the intrinsic shape and detailed dynamics of the galaxy. So 



far, there is little evidence for radial gradients of M/L up 
to 2_Rcff (de Zeeuw 1992, Bertin et al. 1994, CaroUo et al. 
1995, Rix et al. 1997). 

Therefore, we suppose that the galaxy behaves like a 
one-component system. This means that either (I) the stellar 
component is dominant (the invisible matter is negligible), 
or (II) the luminous and dark matters are distributed in 
the same manner, or else (III) the dark matter has a very 
extended distribution, being flat and less important than the 
visible component in the central region (it is the case if the 
luminous part of an elliptical galaxy is embedded in a very 
large dark matter halo). 

The velocity dispersion tensor is not a priori isotropic 
in a galaxy, so that the radial and transverse velocity disper- 
sions may not always be equal (for a discussion see Binney 
& Tremaine 1987). Moreover, normal elliptical galaxies are 
known to have little rotation (cf. de Zeeuw & Franx 1991 
and references therein). We suppose therefore, when using 
the equations of state (^, that the gas of stars is actually 
isotropic and with null angular momentum. This simplifi- 
cation will enable us to calculate the pressure of the gas of 
stars without adding any other parameter to the problem. 

The next step for really computing the entropy is to 
determine the radial mass distribution of galaxies. 



3 THE SERSIC LAW 
3.1 Generalities 

A generalization of the de Vaucouleurs profile has been pro- 
posed by Sersic (1968), Eq. (|^). His profile has been widely 
used in recent years for modelling and theoretical studies of 
elliptical galaxies (e.g., Ciotti 1991; Caon et al. 1993; CoUess 
& Dunn 1996; Courteau et al. 1996; Prugniel & Simien 1997; 
Graham et al. 1996; Graham & CoUess 1997). The Sersic law 
is one of the most compact and simple expressions that de- 
scribes very precisely the light distribution in a wide range 
of early-type galaxies. 

Instead of the historical form of the Sersic law, Eq. (^, 
we will use throughout this paper the following mathemati- 
cally more compact and natural form: 



E(i?) = Eoexp(-(i?/a)") 



(7) 



Sometimes we call it the !/-model, by analogy to the /3-model 
(very often used to describe the x-ray profiles of galaxy clus- 
ters). 

The Sersic profile is characterised by a normalisation 
parameter, Eq; a scale length, a; and a structural (or shape) 
parameter, 1/(1/ = 1/n). The de Vaucouleurs profile is of 
course equivalent to 1/ — 0.25. We call v, a, and Eq the 
primary parameters of the Sersic law. 

From the mathematical point of view, the Sersic profile 
primary parameters are naturally independent, i.e. if - as we 
will show later - after a fit, some correlation among these 
parameters occurs, it will be only due to physical proper- 
ties. This is not the case when other kinds of parameters 
such as effective radius or effective brightness are used, be- 
cause they involve combinations (usually non linear) of the 
primary parameters (cf. Gerbal et al. 1997 and the end of 

Such parameters we call the secondary parameters of 
the Sersic law. 
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We will now give some quantities of astrophysical inter- 
est and/or useful for the present paper. Calculations are gen- 
erally straightforward or are presented in the Appendix 2. 
For the sake of simplicity, we will often give approximate 
but analytical expressions for these quantities. The analyt- 
ical approximations are obtained by fitting convenient sim- 
ple functions to the exact values computed numerically. We 
will refer to the accuracy of the analytical approximations 
of a given quantity C, in terms of (Cthco - Capp)/Ctheo, where 
Cthco is the exact theoretical value and ("app is the analytical 
approximation. Notice that we do not use the calculations 
provided by Ciotti (1991), since he based his calculations on 
a numerical deprojection of the Sersic law while we use an 
analytical approximation described below. 

The total luminosity inside the radius R is readily ob- 
tained by integrating the profile given in Eq. nm: 



n = 0.6213 - 1.25 \og{v) + 0.07656i^ , 



(15) 



L(R) 



(8) 



where 7(6, a:) is the standard incomplete gamma function. 
The total luminosity is, then. 



a Eo -L* (i^) ■ 



(9) 



where L*{u) = 2-kV(2Iv) jv and r(a;) is the complete gamma 
function. 

The total magnitude, mtot = — 2.51ogl/tot+ constant, 
can be obtained accurately (better than 0.5%) for 0.1 <v< 
2.0 by the following analytical approximation: 

-1.44 



m (y) = -0.304;/ - 1.708^" 



(10) 



with m*{u) = — 2.5 log L*(i^). 

The effective brightness, /ieff, is the brightness at the 
effective radius, R^g (defined as the radius having half the 
total luminosity inside it), i.e. ficS = ^i(-Reff), and is given 
by (within 0.5% accuracy): 



Mcff 



-2.5 log Eo + 1.0857(i?cff /a)" ■ 



(11) 



The mean surface brightness inside _Roff can be written, by 
definition, as: 

(^t)eff = m* -^- 5 log(i?eff /a) - 2.5 log Eo + 2.5 log(27r) . (12) 

The effective radius, Rett, is obtained by solving nu- 
merically the equation 2^(2/1/, (R/a)") = r{2/iy). We have 
derived an analytical approximation given by: 

0.6950 - ln(//) 



ln(iicff/a) = 



0.1789, 



(13) 



which is accurate to better than 0.6% in the range 0.1 < 
u < 1.8 for f?efr/a- 

Another useful quantity is the Petrosian metric radius, 
which is related to the logarithm slope of the brightness 
profile (Petrosian 1976). The Petrosian radius is defined with 
the help of the variable 



V -- 



-2.5 log 



dlog L{R) \ 
2d log i? J 



M(i?) - MR)) 



(14) 



where /i(-R) is the magnitude (equal to —2.5 log(E)) at R and 
{fj.{R)) is the mean magnitude inside R (Sandage & Perel- 
muter 1990). For a given value of rj, a radius Rp is defined. 
Besides being model independent, the Petrosian radius gives 
a metric measure that is distance independent. For the par- 
ticular case where Rp = ReS, the Petrosian parameter rj is 
a function of v given approximately by: 



which is valid in the range 0.1 < v < 1.8 with an accuracy 
of 0.1%. 

As we mentioned before, almost all 'astrophysical' quan- 
tities such as the effective radius, the effective brightness and 
the mean effective brightness are non-linear combinations of 
the 'primary parameters' Eo, a and v. 

It is also interesting to rewrite the Sersic profile (Eq. 0) 
introducing explicitly the effective Radius: 



E(i?) = Eoexp[-64i?/i?,ff)''] 



(16) 



with bi, = [R^ft/aY ■ The original Sersic profile, Eq. (|l|), is 
obtained making fc„ = 1.08576^ and iJt{R) = — 2.5 log E(_R). 



3.2 Three dimensional Sersic deprojection 

The two dimensional Sersic density distribution has no an- 
alytical deprojection in three dimensions. For our purposes, 
numerical deprojection (cf. Young 1976 for the de Vau- 
couleurs case) or the asymptotic behaviour (Ciotti 1991) 
of the Sersic profile are not well adapted. Mellier & Mathez 
(1987) have given a simple analytical approximation of the 
deprojection of the de Vaucouleurs profile; we will use a 
modified form of the Mellier & Mathez three dimensional 
density profile, generalised for the Sersic profile case: 



p{r) = po 



oxp 



with the normalisation: 



Po = Eo 



r(2/:.) 



M 



2ar((3-p)/i/) L 



(17) 



(18) 



which is obtained by requiring that both total masses cor- 
responding to Eqs. (§) and ^ be equal. M/L is the mass 
to luminosity ratio and we suppose hereafter that it is inde- 
pendent of t he g alaxy radius. This is the above mentioned 
hypothesis (§2.2) that the galax;y behaves like a one compo- 
nent system. Notice, however, that up to this point we do 
not require that all galaxies have the same M/L ratio. 

The exponent p in Eq. (|l^) is obtained by fitting the 
(numerical) Abel deprojection of the two dimensional profile 
Eq. (Q) to the analytical three dimensional profile Eq. ( p^ . A 
high quality match, better than 5%, between the numerical 
deprojection of Eq. (^ and Eq. ( [l7| ) is obtained when: 

p = 1.0 - 0.60971/ + 0.05463i/^ , (19) 

for 0.1 < V < 1.8, and the radial range 10~^ < R/Rctt < 
10''. The discrepancy between the numerical and analytical 
deprojections is more pronounced for very small and very 
large radii. Prugniel & Simien (1997) proceeded in the same 
fashion, using a modification of the Mellier & Mathez profile, 
but their deprojection is only valid in the range 0.1 < u < 
0.5. Within this range their result and ours agree well within 
1%. For 0.5 < < 1 our result is within 5% to Prugniel & 
Simien's, while for v > 1, our result diverges parabolically 
from theirs. 

Notice that the behaviour of Eq. ( p^ as r — » is not 
the same as the exact Abel deprojection of the Sersic profile. 
For small radii, pAbei{r) oc r""^ for v < 1. This should not 
introduce any appreciable error in our analytical approxi- 
mation; the fraction of mass inside 10~^ < R/Rc« < 10"' is 
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97% for V = 0.1, 99.6% for u = 0.25, and greater than 99.9% 
for u > 0.55. 

The corresponding mass inside the radius r, M(r), is 
readily obtained integrating Eq. (|l7|) in a spherical volume: 



M{r) = — — 7 ' ^ 



(20) 



Additionally, we give the relation between the two di- 
mensional half-mass radius R^s and the three dimensional 
half- mass radius, reft, as a polynomial approximation: 

r^s/R^s = 1.356 - 0.0293u + 0.0023;/^ , (21) 

This approximation agrees well within 0.25% to the values 
computed by Ciotti (1991) who used the numerical Abel in- 
version of the Sersic profile. This reassures us on the quality 
ofEq. (|l]). 

We define a gravitational radius by the relation rg = 
GM^ /\U\, with U the gravitational potential energy. Again, 
we give a simple analytical approximation: 

0.820 - 0.924 ln{u) 



ln(rg/a) = 0.845 + 



(22) 



that is accurate to better than 0.8% in the range 0.1 < v < 
1.8 for Vg/a. This radius is a combination of the natural scale 
length and the form factor; this is naturally due to the loss 
of homology when moving from a de Vaucouleurs law to a 
v-mode\. A similar behaviour occurs with -Reff, Eq. (^. It 
turns out that the relation between _Rcff and Vg is roughly 
linear, independent of given approximately by 



ln(rg/a) = 1.16 + 0.98 ln(i?cff /a) : 



(23) 



accurate to better than 5% in the range 0.1 < v < 1.8 for 
rg/a. For 0.15 < < 1.0, the relation is practically linear 
(within 15%), viz. rg ~ 3.0i?cff ~ an interesting result given 
that the Sersic model is not homologous. 



3.3 Entropy and the Sersic profile 

Given a mass distribution law, the deprojected Sersic law, 
we can make the connection between the specific entropy 
and the parameters of the light profile. We need, however, 
an additional hypothesis to get the pressure corresponding 
to Eq. (^^. We consider that the pressure is isotropic and 
the galaxy is in hydrostatic equilibrium. This implies an 
isotropic velocity dispersion tensor and the simple relation, 
P{r) = p{r)a^, relating the pressure, density, and radial ve- 
locity dispersion, ar- With the isotropy hypothesis, ar (or 
P{r)) is, thus, obtained as a solution of the hydrostatic equi- 
librium equation: 



— [p{r) a^rir)] 



G M{r) p{r) 



(24) 



With the above hypothesis, one may substitute the den- 
sity and pressure in Eq. (^ and a relation of the kind 
s{a, Eq, v) can be derived. It turns out that s{a, Eq, u) may 
be well approximated by the analytical expression: 

1 5 
s(a, Eo,!/) = -ln(Eo) + -ln(a) 

+ ci ln(i/) - i +02!^"=* + CO . (25) 

The constants Ci as well as the derivation of Eq. (^) are 
given in Appendix 2. 



It is interesting to note that for a Plummer profile (poly- 
trope of index 5, e.g. Binney & Tremaine 1987), where we 
are able to compute exactly the specific entropy, we have 
the same dependency of s on the Plummer scale and nor- 
malisation parameters as in Eq. (p5[). 

The relations that can be written as the formal expres- 



s(a, Eo, = Si 



(26) 



where Si means any given value of the specific entropy, define 
a family of surfaces in the three dimensional space defined by 
the parameters Eq, a, v. By analogy with the Fundamental 
Plane, we call such a surface the 'Entropic Plane'. 

We do not know a •priori the value of the specific en- 
tropy. We can, nevertheless, propose two hypotheses con- 
cerning the entropy of elliptical galaxies: 

• Hypothesis I: The specific entropy of early- type galaxies 
in a given cluster is constant. This is the weak hypothesis. 

• Hypothesis II: The specific entropy of early-type galax- 
ies in clusters is universal, independent of the cluster in 
which they are located. This is the strong hypothesis. 

One may hope that in the future, some theory of galaxy 
formation will move one of the two hypotheses above into 
a firm knowledge and provides us with the true value of s. 
In the absence of a way to calculate s from first principles, 
these hypotheses must be tested using the observed bright- 
ness profiles of early type galaxies. 



4 THE DATA 

In order to verify these hypotheses, we use the brightness 
profile data of elliptical galaxies belonging to two rich clus- 
ters of galaxies - Coma and ABCG 85 - and the group 
associated with NGC 4839 (which is in the neighbourhood 
of Coma). 



4.1 Description 

Galaxies in the Coma cluster have been selected from the 
catalogue of Biviano et al. (1995a) among those in the cen- 
tre of the cluster and in the group around NGC 4839, being 
brighter than 17.8 V mag, fainter than 13 V mag (to exclude 
the biggest, saturated galaxies), with elliptical morphologi- 
cal types as assigned by Godwin, Metcalfe & Peach (1983, 
hereafter GMP), and with measured redshifts that confirm 
them as cluster members (radial velocity in the range 3000- 
10000 km s"\ cf. Biviano et al. (1995b)). The elliptical 
galaxies of ABCG 85 were extracted from the catalogue of 
Slezak et al. (1998), by selecting those galaxies belonging to 
the cluster (from available redshifts, with velocities in the 
range 13350-20000 km s"\ cf. Durret et al. (1998)) and 
with elliptical appearance in the CCD images. The final list 
of elliptical galaxies used for the fitting is given in Tables |l|, 
^ and ^ (in their first columns we give the corresponding 
identifications in both catalogues, Biviano et al. 1995a for 
Coma and the NGC 4839 group, and Durret et al. 1998 for 
ABCG 85). 

We have used V images that had been reduced and cal- 
ibrated following standard techniques (see Lobo et al. 1997, 



© 1999 RAS, MNRAS 000, 



6 G.B. Lima Neto et al. 



Slezak et al. 1998). We then cleaned the images by eliminat- 
ing foreground stars that contaminated our target galaxies. 
This was made through specific tasks within IDL (Interac- 
tive Data Language, Research Systems, Inc.). The elimina- 
tion was done by interpolating the field around the star by a 
4th degree polynomial. We then obtained the growth curves 
(integrated magnitudes within circular regions) by means of 
the 'ellipse' task in 'stsdas' environment of irafQ For every 
galaxy in our sample we obtained a table listing the radii 
of the circular apertures and the corresponding integrated 
magnitudes inside the apertures. Galaxies that had less than 
nine data points were rejected. This means that the galax- 
ies on the faint surface brightness end were rejected by our 
selection. Since our analysis is not sensitive to the complete- 
ness of the sample, the rejection of very faint galaxies should 
have no effect in our results. 

The use of a circular aperture was motivated by two rea- 
sons. First, in order to be able to compare directly our results 
with those found in the literature, in particular Graham & 
CoUess (1997) and, second, because the theory we developed 
is based on spherical symmetry of the mass distribution. No- 
tice that previously (Gerbal et al. 1997) we have employed a 
different method for obtaining the growth curves, taking into 
account the shape of the isophotes (we computed the total 
luminosity inside a given isophote) . Comparing our previous 
results (Figs. 1 and 2 from Gerbal et al. 1997) with those 
presented here shows that the difi'erence is not significant 
(i.e., the pairwise distributions of the Sersic parameters do 
not change). 



4.2 The fitting method 

We have fitted the resulting luminosity growth curves of 
the elliptical galaxies in Coma (including the group of 
NGC 4839) and ABCG 85. 

The advantage of fitting the integrated fiux instead of 
the luminosity profile is that the latter is more sensitive to 
irregularities in the light distribution. Besides, fitting the 
light growth curve gives more weight to the inner than the 
outer parts of the galaxy, refiecting better the galaxy overall 
structure (D'Onofrio & Prugniel 1997). Therefore, since the 
inner region of a given galaxy is probably more relaxed than 
the outer region (which is more susceptible to perturbations) 
it seems more reasonable to give more weight to the inner 
region when computing the specific entropy. The drawback is 
that if the galaxy has some unusual centre, an active nucleus, 
star burst or some exceptional darkening, then these features 
will be propagated in the whole growth curve. In this case, 
fitting the luminosity profile would be better to subtract 
the anomalous central contribution (D'Onofrio & Prugniel 
1997). Moreover, Burstein et al. (1987) pointed out that the 
fitting of growth curve have systematic uncertainties that 
may lead to correlated errors of the total magnitude and the 
effective radius. This last point may be minimized by fitting 



* IRAF is the Image Analysis and Reduction Facility made avail- 
able to the astronomical community by the National Optical As- 
tronomy Observatories, which are operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., under 
contract with the U.S. National Science Foundation. 




V 



Figure 1. Distribution of the shape parameter u for all galaxies, 
grouped by cluster membership. The distribution does not peak 
at 0.25, the value corresponding to the de Vaucouleurs profile. 

t he uncorrelated primary parameters of the Sersic law (cf. 
^.l[ ) rather than their non linear combinations. 

Consequently, we use the integrated form of the bright- 
ness profile as given by Eq. (^). The luminosity growth 
curve fits are done with a standard least square minimisa- 
tion method, using the 'minuit' programme from the CERN 
software library. In order to reduce the effects of the seeing 
(FWHM « 0.9 arcsec for Coma and 1.2 arcsec for ABCG 85) 
and/or eventual nuclear activity, the data points are taken 
only from 2.5 arcsec outwards. For galaxies of the same clus- 
ter, the effect of using a central exclusion mask fixed in an- 
gular units, rather than a linear measure, is virtually negli- 
gible (in this precise case the galaxies can be considered to 
be at the same distance from us; for Coma the error should 
be smaller than ~ 5%). A difficulty may rise in comparing 
galaxies in different clusters. However, in our case, the clus- 
ters are not more than a factor 2.5 apart (radial distance; cf. 

below) and changing the size of the mask from 0.6 to 3.0 
arcsec did not alter significantly the fitting results, i.e., the 
values obtained for the Sersic parameters were compatible 
within their 3(t error bars [notice that this test on the size 
of the central mask were also done by Gerbal et al. (1997)]. 

We have used the growth curve up to the surface bright- 
ness isophote at 24.0 V mag/arcsec^, corresponding to a 3o- 
signal above the background sky. This is done in order to 
minimize errors introduced in the background subtraction. 
We performed the fits varying the outer cut-off from 22.5 to 
24.5 V mag/arcsec^; the fitting results were still basically 
the same, i.e., the values obtained for the Sersic parameters 
were compatible within their 3(t error bars. 



5 FITTING RESULTS AND CORRELATIONS 

The results of the growth curve fits are shown in Tables |l|, ^ 
and I for Coma, ABCG 85, and the group of NGC 4839, re- 
spectively. We present also in these tables the usually quoted 
effective radius, RcS, and mean surface brightness {/i)cff, 
both computed with formulas ( |l^ ) and (|l2|), respectively. 
The errors in the Sersic parameters are represented graphi- 
cally in the Figures |^-^. 

In Fig. n we show the distribution of the shape param- 
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Table 1. Growth curve fitting results for Coma. Ident is the GMP 
number and name is either an IC or NGC identification, v and 
a are the shape and scale parameters fitted with the Sersic law. 
/iQ = — 2.51og{So) is the intensity parameter and the magnitude 
is in the V band. iJcH and i a\„a are calculated from the primary 
parameters as explained in § p.l| 

Ident* name v log(a) /io -Rcff (A')eff 

(arcsec) (mag/D") (arcsec) (mag/D") 



* From the catalogue of Biviano et al. (1995a) 
© 1999 RAS, MNRAS 000, 0-|l| 
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Table El Continued 



Ident* 


name 


u 


log (a) 


MO 




(M>cff 








(arcsec) (: 


mag/n") 


(arcsec) 


(mag/n") 


3656 




0.20 


-3.55 


13.65 


16.07 


22.39 


3664 


NGC4864 


0.49 


-0.55 


16.25 


4.26 


19.28 


3681 




0.62 


-0.34 


18.92 


2.53 


21.13 


3707 




0.47 


-0.74 


18.02 


3.31 


21.19 


3733 


IC3960 


0.25 


-2.79 


13.27 


5.67 


20.21 


3782 




0.33 


-1.79 


15.05 


3.38 


20.06 


3792 


NGC4860 


0.45 


-0.74 


15.58 


4.03 


18.91 


3794 




0.50 


-0.87 


16.59 


1.89 


19.56 


3855 




0.63 


-0.20 


19.54 


3.36 


21.71 


3914 




0.56 


-0.70 


16.11 


1.68 


18.67 


4103 




0.74 


-0.02 


19.43 


3.07 


21.16 


4129 




1.04 


0.07 


19.88 


1.86 


20.94 


4200 




0.55 


-0.61 


16.87 


2.15 


19.47 


4230 




0.23 


-3.05 


12.81 


7.90 


20.36 


* From the catalo 


guo of Biviano et al. (1995a) 






Table 2. Growth 


curve 


fitting results for ABCG 85 The columns 


have the same meaning as in Table |l| 


, except for the Ident column. 


Ident* 


u log (a) 


MO 


^off 


(M)cff 






(arcsec) 


(,mag/u ) 


(arcsec) ( 


mag/ u ) 




152 


0.59 


-0.63 


17.9 


1.56 


20.30 




156 


0.20 


-4.40 


13.0 


3.52 


22.00 




175 


0.29 


-2.49 


14.7 


2.14 


20.56 




179 


0.50 


-0.92 


17.7 


1.58 


20.68 




182 


0.69 


-0.02 


18.6 


3.69 


20.50 




197 


0.38 


-1.31 


15.8 


3.05 


19.95 




202 


0.71 


-0.30 


18.2 


1.81 


20.05 




208 


0.26 


-2.61 


15.5 


5.91 


22.20 




209 


0.76 


-0.07 


18.9 


2.53 


20.60 




212 


0.64 


-0.58 


18.3 


1.30 


20.42 




214 


0.46 


-0.62 


18.1 


4.94 


21.40 




218 


1.16 


0.12 


20.6 


1.79 


21.52 




222 


0.63 


-0.54 


18.5 


1.51 


20.67 




225 


0.71 


-0.32 


18.7 


1.74 


20.54 




228 


0.89 


-0.07 


20.2 


1.78 


21.48 




229 


0.80 


-0.21 


20.1 


1.66 


21.71 




231 


0.46 


-1.03 


20.5 


1.82 


23.78 




235 


0.49 


-0.87 


17.7 


1.93 


20.74 




236 


0.44 


-0.96 


16.3 


2.70 


19.77 




246 


1.43 


0.31 


22.8 


2.18 


23.40 




253 


0.62 


-0.48 


18.3 


1.77 


20.49 




263 


0.24 


-3.59 


13.1 


1.49 


20.36 




283 


0.34 


-1.49 


16.8 


5.35 


21.67 




296 


0.85 


0.01 


19.6 


2.30 


20.98 




305 


0.76 


-0.19 


19.4 


1.93 


21.05 




310 


0.94 


-0.12 


19.8 


1.43 


21.04 




316 


1.05 


0.25 


19.9 


2.77 


20.92 




324 


0.35 


-1.31 


17.9 


6.27 


22.54 




329 


0.66 


-0.42 


19.4 


1.70 


21.39 




* From the catalogue of Durret et al. (1998) 







eter v for all our fits. The distribution is fairly regular and 
unimodal, but it is not symmetric (skewness ~ 1). It is in- 
teresting to note that u = 0.25 - the value corresponding to 
the de Vaucouleurs profile - is neither the median {v = 0.50) 
nor the mean {f — 0.55) of the distribution. The maximum 
of the distribution corresponds to ~ 0.4. Very high values 
of u are obtained for a few galaxies. Although the galaxies 
selected by us are classified as ellipticals, our fits suggest 
that some of them may be in fact either lenticulars or dwarf 
spheroidal galaxies (those with u ^ 1). Notice, however, that 



2727 


IC4026 


0.27 


-2.27 


14.33 


6.74 


20.58 


2736 




0.71 


-0.15 


19.31 


2.51 


21.14 


2753 




0.79 


0.00 


19.69 


2.73 


21.27 


2777 




0.53 


-1.01 


16.68 


1.02 


19.41 


2794 


NGC4898B 


0.49 


-0.65 


16.48 


3.29 


19.50 


2798 


NGC4898A 


0.50 


-0.60 


15.84 


3.28 


18.76 


2805 




0.32 


-2.04 


14.27 


2.42 


19.46 


2839 


IC4021 


0.32 


-1.91 


14.14 


2.94 


19.24 


2879 




0.69 


-0.15 


19.38 


2.73 


21.28 


2897 




0.64 


-0.13 


18.84 


3.71 


20.97 


2910 




1.09 


0.23 


18.82 


2.48 


19.79 


2910 




0.50 


-0.62 


16.67 


3.23 


19.61 


2914 




0.43 


-0.67 


18.35 


6.42 


21.92 


2921 


NGC4889 


0.45 


-0.32 


16.10 


10.86 


19.45 


2922 


IC4012 


0.53 


-0.66 


15.96 


2.30 


18.71 


2940 


IC4011 


0.34 


-1.53 


15.37 


4.05 


20.07 


2960 




0.43 


-0.80 


17.53 


4.62 


21.09 


2975 


NGC4886 


0.30 


-1.97 


14.35 


5.48 


19.99 


2976 




0.63 


-0.34 


19.10 


2.34 


21.25 


3058 




1.07 


0.34 


20.46 


3.31 


21.48 


3073 


NGC4883 


0.44 


-0.81 


16.10 


4.19 


19.59 


3084 




0.42 


-1.00 


16.27 


3.24 


19.90 


3113 




0.66 


-0.21 


19.15 


2.73 


21.18 


3126 




0.87 


0.10 


19.56 


2.73 


20.93 


3133 




0.46 


-0.91 


17.19 


2.64 


20.50 


3170 


IC3998 


0.33 


-1.52 


15.52 


5.55 


20.43 


3170 


IC3998 


0.34 


-1.47 


15.50 


4.89 


20.24 


3201 


NGC4876 


0.52 


-0.53 


16.65 


3.36 


19.46 


3206 




0.54 


-0.47 


17.56 


3.19 


20.21 


3213 




0.42 


-1.07 


15.95 


2.98 


19.64 


3222 




0.48 


-1.00 


15.73 


1.61 


18.82 


3222 




0.93 


-0.08 


17.61 


1.59 


18.85 


3269 




0.34 


-1.83 


14.87 


2.50 


19.73 


3291 




0.33 


-1.42 


16.92 


7.26 


21.86 


3292 




0.64 


-0.23 


19.20 


2.97 


21.34 


3296 


NGC4875 


0.39 


-1.29 


15.17 


2.92 


19.22 


3302 




0.45 


-0.93 


17.32 


2.66 


20.66 


3312 




0.51 


-0.89 


18.11 


1.60 


20.99 


3329 


NGC4874 


0.60 


0.20 


17.70 


10.01 


20.01 


3340 




1.07 


0.13 


20.77 


2.04 


21.77 


3352 


NGC4872 


0.50 


-0.72 


15.77 


2.56 


18.71 


3367 


NGC4873 


0.35 


-1.40 


15.40 


5.19 


20.07 


3367 


NGC4873 


0.30 


-1.82 


14.80 


6.67 


20.32 


3400 


IC3973 


0.34 


-1.82 


13.65 


2.55 


18.51 


3414 


NGC4871 


0.23 


-3.17 


12.63 


6.12 


20.19 


3423 


IC3976 


0.32 


-2.01 


13.61 


2.59 


18.78 


3439 




0.62 


-0.19 


18.60 


3.56 


20.81 


3486 




0.78 


-0.26 


18.79 


1.55 


20.42 


3487 




0.41 


-1.27 


15.94 


2.29 


19.78 


3489 




0.91 


-0.04 


19.99 


1.83 


21.29 


3510 


NGC4869 


0.47 


-0.66 


16.02 


4.24 


19.25 


3522 




0.43 


-1.14 


15.61 


2.14 


19.17 


3534 




0.58 


-0.43 


18.23 


2.71 


20.67 


3554 




0.37 


-1.33 


16.80 


3.99 


21.16 


3557 




0.51 


-0.48 


17.54 


3.98 


20.40 


3561 


NGC4865 


0.43 


-0.89 


15.15 


3.70 


18.69 


3565 




0.65 


-0.23 


19.85 


2.84 


21.94 


3639 


NGC4867 


0.54 


-0.57 


16.03 


2.59 


18.70 
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Table 3. Growth curve fitting results for the N4839 group. The 
columns have the same meaning as in Table ^ 



Ident* 


name 


V 


log(a) 
(arcsec) 


(mag/D") 


_Reff 

(arcsec) 


(M>cff 

(mag/n") 


4714 




0.48 


-0.83 


17.29 


2.38 


20.38 


4792 


NGC4842B 


0.52 


-0.62 


16.65 


2.62 


19.42 


4794 


NGC4842A 


0.49 


-0.61 


16.05 


3.75 


19.10 


4852 




0.69 


-0.32 


18.60 


1.84 


20.49 


4907 




0.41 


-1.17 


15.26 


2.75 


19.06 


4928 


NGC4839 


0.25 


-2.03 


14.63 


27.87 


21.46 


4937 




1.07 


0.26 


20.60 


2.78 


21.62 


5051 




0.32 


-1.77 


14.17 


3.89 


19.24 


5102 




0.57 


-0.22 


19.07 


4.65 


21.56 



* From the catalogue of Biviano et al. (1995a) 




10-2 
a (arcsec) 



Figure 2. Distribution of the intensity parameter —2.5 log So (in 
V magnitude) and the scale length a (in arcsec) for all galaxies, 
each panel for a cluster. The bars arc \u errors. 



our sample is not complete (we had to reject a few faint 
galaxies, the brightest saturated ones, and we selected only 
the ones with known redshift) and the distribution will be 
underestimated for the low surface brightness galaxies (i.e. 
high v\ 

In Figs. ^ ^, and ^ we display the relations between the 
primary parameters So, a, and i/. It is obvious that all these 
quantities are strongly pairwise correlated, although non- 
linearly. Before analysing quantitatively these correlations, 
let's make the connection with the Entropic Plane. 

From Eq. (^^, we expect that there are indeed non- 
linear correlations among the Sersic profile parameters, since 
that equation defines a surface in a three dimensional space. 
In fact non-linear correlations are expected because the En- 
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Figure 3. Distribution of — 2.51ogSo(in V magnitude) and v. 
The bars are Icr errors. 



tropic Plane relates non-linearly these parameters (if the 
weak or strong hypotheses hold) . In order to verify whether 
or not Eq. (|2^) actually holds, we have defined two auxiliary 
variables: 



y(So,a) = -ln(Eo) 



and 



CO + ci \n{v) - - 



(27) 



(28) 



(the constants d are defined in Appendix 2). Notice that 
X{y) is a monotonic descending function of v only. With 
the above definitions, we can rewrite the specific entropy as 
a simple linear relation as: 



y(So,a) 



So 



X{u). 



(29) 



This relation defines a family of straight lines - 'En- 
tropic Lines' which correspond to the Entropic Plane seen 
edge-on - for which the slope is —1 and the zero-point value 
(i.e. So, the specific entropy) is unknown. 

It is straightforward to calculate Yi and Xi from the val- 
ues of So, a, and u for the galaxy i and, consequently, the 
specific entropy, so for each galaxy (but notice that so cal- 
culated in this way is distance dependent - we will develop 
this in §^). In Table ^ we give for each cluster the mean 
value of So, as well as the standard error and deviation. The 
zero-point of the Entropic Lines is simply the mean value of 
So shown in Table ^. 

In Fig. P we display Yi and Xi for all our galaxies; the 
Entropic Line is superimposed to the data points, where we 
have used the values quoted in Table ^. We also show the 
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1 1.2 1.4 



Figure 4. Distribution of a (in arcsec) and u. Notice that liere 
we have the best correlation between two of the three parameters 
of the Sersic profile. The bars are Icr errors. 



Table 4. Statistic of the specific entropy. The mean value of the 
specific entropy, sqi is used as the zero-point of the Entropic Line 
(cf. Fig. ^). While the Standard Deviation refiects the individual 
scatter, the Standard Error takes into account the number of data 
points. 





Coma 


Abcll85 


NGC4839 


Mean so 


-4.63 


-5.99 


-4.28 


Standard Error 


0.14 


0.20 


0.61 


Standard Deviation 


1.20 


1.09 


1.84 


Number of galaxies 


73 


29 


9 



displacement vector for a typical error in v (hvjv = 0.05); 
this is done taken in account the correlations of v with both 
So and a (cf. Figs ^ and ^, as well de definition of 

The galaxies follow remarkably the predicted linear cor- 
relation, in spite of all our simplifying hypotheses. The small 
dispersion around the line representing Eq. (|2^) strongly 
suggests that in each individual cluster, these galaxies share 
the same specific entropy: at least the weak hypothesis is 
fulfilled. In Fig. |, we show the residual distribution, i.e., 
[Yi — yi,thcor)/yi,thoor (whcrc Yi IS obtained from the obser- 
vational data and we use Eq. ( psj ) to obtain yi,theor from 
So — Xi); the standard deviation is equal to 9.5%. Applying 
the Kolmogorov-Smirnov test to the residual distribution 
shown in Fig. ^ results that its distribution is compatible to 
a gaussian one. Notice, however, that due to the small num- 
ber of data points, only a very deviant distribution would 
not be compatible with a gaussian. 
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ABCG 85 




++ + ++ 



NGC 4839 



5 10 15 20 25 

X{v) 

Figure 5. The specific entropy relation and their corresponding 
residual for the elliptical galaxies in Coma, ABCG 85, and the 
group of NGC 4839. The straight lines are the 'Entropic Lines', 
(Eq. (pgj), where the slope is equal to —1 and the zero-points, sq, 
are from Table |: -5.99 ±0.20, -4.63 ±0.14, and -4.28 ±0.61 for 
ABCG 85, Coma and NGC 4893, respectively. Note that X{v) is 
a monotonic descending function of i.e., high values of X{v) 
correspond to small u. The residual plot below each Entropic 
Line is ^data ~ ^predicted- The arrow in the Coma panel shows 
the typical displacement vector for an error in u of 5%. This 
displacement is due to the definition of X{v) and the correlations 
with So and a. 



It would be possible to reduce the scatter of Xi and Yi 
by using a linear fit (slope and zero-point as free parameters) 
or a polynomial one instead of computing so the way we did 
and using the predicted Entropic Line. This, however, would 
bury the physics behind the Entropic Plane relation - in fact, 
the variables X(v) and y(Eo,a) would lose their meaning. 
Anyway, fitting the Xi and Yi data points to a free-slope 
straight line (using a standard least-square method) yields 
the following slopes: -0.84±0.05, -0.93±0.06, and -0.66± 
0.12 for Coma, AbeU 85, and NGC 4839 group, respectively. 
The mean value of the slope, taking into account the error 
bars, is 0.86 ± 0.14. Although this is statistically consistent 
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NGC 4839 
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Figure 6. Dispersion distribution of the galaxies around the 'En- 
tropy Line' shown in Fig. ^. This error is associated with the un- 
certainty in determining the specific entropy of each galaxy (see 
§^ and the relative distances between clusters (see ^). We show 
also the best fit (standard least-square method) gaussian to all 
data points and the one standard error, cr, in percent. 
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Figure 7. Scatter diagram of v and -Reft- Notice the weaker cor- 
relation in these plots compared to the correlations between the 
primary parameters v and a. 



with the predicted Entropic Line, the tilt is almost certainly 
real. There is a number of effects that may account for the 
scatter and the deviation of the data points from a hne of 
slope exactly equal to —1. This point will be addressed in 

^ n 

Similar to our Fig. M, Young & Currie (1995) have al- 
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Figure 8. Scatter diagram of {fJ.)eS and -Refl- The dotted lines 
are from the relation logij^ff 0.32(^)efj found by Oegerle & 
Hoessel (1991) for brightest cluster members. Notice the weaker 
correlation in these plots compared to the correlations between 
the primary parameters Sq and a. 



Table 5. Non parametric correlation tests between the param- 
eters of the Sersic-law for the Coma cluster. Both p and r are 
defined in the interval [—1,1], the value meaning no correla- 
tion. Higher absolute values of z indicate a greater significance of 
the correlations. 



Parameters 


IpI 


Spearman 


2 


r Kendall 


2 


X(i.)-Y(Eo,a) 




0.94 


7.97 


0.82 


10.3 


u-a 




0.95 


8.10 


0.84 


10.5 






0.89 


7.57 


0.73 


9.18 


a-Eo 




0.89 


7.56 


0.73 


9.18 






0.51 


4.36 


0.37 


4.65 






0.18 


1.54 


0.13 


1.60 


Table 6. Same as 


Table || for ABCG 85. 






Parameters 


IpI 


Spearman 


2 


r Kendall 


2 


X{u)-Y{T:o,a) 




0.95 


5.03 


0.84 


6.38 


v-a 




0.97 


5.13 


0.88 


6.68 






0.89 


4.71 


0.80 


6.04 


a-Eo 




0.85 


4.52 


0.74 


5.63 






0.32 


1.67 


0.22 


1.65 






0.31 


1.64 


0.21 


1.61 
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Table 7. Same as Table || for NGC 4839 group. 



Parameters 


\p\ Spearman 


z 


|r Kendall 


2 


X{iy)-Y{SQ,a) 


0.97 


2.73 


0.89 


3.34 


v-a 


0.97 


2.73 


0.89 


3.34 




0.92 


2.59 


0.78 


2.92 


a-So 


0.92 


2.59 


0.78 


2.92 




0.42 


1.18 


0.33 


1.25 




0.20 


0.56 


0.06 


0.21 



ready presented the v-a correlation for dwarf ellipticals in 
the direction of the Fornax cluster (figure 1 in their paper) 
and fitted it with a polynomial having five free parameters. 
Their small r.m.s. scatter (0.108) suggested that the u-a cor- 
relation could be a useful distance indicator. More recently, 
Binggeli & Jerjen (1998), who have also fitted the growth 
curve using the Sersic profile, have presented the same cor- 
relation between the shape parameter and the scale length 
for early-type galaxies belonging to the Virgo cluster (Fig. 7 
in their paper). They have fitted a parabola to the v-a cor- 
relation and have derived a r.m.s. equa l to 17%. 

We have mentioned in Section 3.1 that the 'astrophysi- 



cal' parameters are non-linear combinations of the primary 
parameters, ft is therefore not surprising that the correla- 
tions between these quantities are less clear compared to 
the corresponding relations obtained when one uses the pri- 
mary parameters. For instance, the distribution of 7?eff and 
V (Fig. 1^ compared to the distribution of a and v (Fig. ^ 
shows how the primary parameters correlates better than 
the classical astrophysical parameters. One can also see the 
Fig. 3 from the paper by Graham et al. (1996), where they 
plot the ReS-v relation and compare it with the a-u rela- 
tion displayed in Fig. 7 from the paper of Binggeli & Jerjen 
(1998) for Virgo galaxies. Finally, one can compare the com- 
puted values of -Reff and (^)off (Fig. ^) and compare it to 
the relation between —2.5 log So and a (Fig. |^. 

As already noted by Oegerle & Hoessel (1991), the scat- 
ter is rather larger for normal elliptical galaxies than for the 
brightest cluster members (BCG). fn Fig. ^ we have super- 
posed the empirical relation logi?eff tx 0.32(/i)efr found by 
Oegerle & Hoessel (1991, their Fig. 5) for their sample of 
BCGs. It is clear from Fig. ^ that this is the most scattered 
relation. 

The relations observed in Figs. |^, ^, and ^ are nothing 
more than projections of the Entropic Plane predicted by 
Eq. disl ) onto the three planes defined by the parameters 
Eo, a, and u. The thickness of the pairwise correlations seen 
in these figures is largely due to the fact that the specific En- 
tropic Plane is not seen exactly edge-on. The non-linear com- 
binations of these parameters (the ReB-{ij.)cB correlation, for 
instance) present, as one would expect, greater scatter. By 
the same token, one could naively expect that the scatter of 
the Xi-Yi relation (which is also a non-linear combination 
of the Sersic parameters) would be greater than the scatter 
of the primary parameters pairwise correlations. 

In Tables ^ to ^, we quantify the quality of the Sersic pa- 
rameters pairwise relations using two non-parametric anal- 
ysis: the Spearman and the Kendall rank correlation tests 
(for a detailed description, see Siegel & Castellan 1989; Press 
et al. 1992). (In ^ below, we will give the scatter around 
the Entropic Line in terms of the distance modulus for each 



galaxy in Fig. ^) As it can be seen, the relations involv- 
ing the secondary parameters, RcH and {/i)cff , have indeed 
smaller correlation coefficients than the relations involving 
only the primary Sersic parameters. 

The Entropic Line turns out to have about the same 
correlation coefficients as the iz-a relation (the latter being 
slightly more tightly correlated). This tells us that the En- 
tropic Line does not seem to be just an arbitrary non-linear 
combination of the primary parameters, but may be indeed 
a fundamental relation based on the assumption of unique 
specific entropy. It also tells us that the i/-a relation is al- 
most an edge-on projection of the Entropic Plane. The fact 
that the Entropic Line has not a better correlation coefficient 
may be due to one or some of our simplifying hypotheses in 
deriving the specific entropy relation for elliptical galaxies 
(for instance, the ideal gas, the isotropic pressure, and the 
mass-luminosity ratio, see also gTI). 



6 TOWARDS A DISTANCE INDICATOR 
6.1 Distance dependence 

The specific entropy can be written as a function of the 
three parameters of the Sersic profile (cf. Eq. If we could 
somehow determine the value of the specific entropy for a 
given galaxy, then we would have a relation linking one of 
the parameters as a function of the remaining two others, 
for instance, u(a,T,o). 

Any relation that operates between a dimensionless pa- 
rameter and one or more distance-dependent parameters can 
be potentially used as a distance indicator. In the present 
case, a correlation between the distance independent v and 
the distance dependent a and Eq could be used as a yard- 
stick. 

However, a difficulty arises since we need to know a 
priori the value of the specific entropy. This problem can 
be solved partially if we consider relative distances between 
clusters of galaxies (but not individual galaxies) instead of 
absolute distances. Measuring the distances of clusters in- 
stead of individual galaxies reduces the error induced by the 
dispersion in the specific entropy relation. 

The parameters a and Eq are, so far, given in metric 
units, not angular. The distance dependence of Y, defined 
in Eq. (p7[), can be written explicitly by converting a and 
Eo to those observed, a" and Eq", as follows: 



y = -hj:o") 



^ln(a") + ^ln(i3) 



(30) 



where D is the distance to the cluster. 

The fits of the growth curves provide us with the pa- 
rameters Eo", a", and u, thus the parameters X and Y 
(without the distance term). Given a sample of early type 
galaxies from a cluster of known distance, we obtain the 
specific entropy as the mean value of so found for all sample 
galaxies in that cluster, provided that Hypothesis I holds. 



6.2 Relative distances 

Likewise, given two samples of galaxies on two different clus- 
ters of unknown distances, and assuming that Hypothesis II 
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holds, we could determine the relative distance of the clus- 
ters directly from Eq. (]30|). The relative distance between 
two clusters is given by 



(31) 



where Si's are the mean entropy of each cluster. To these 
distance corresponds the relative distance modulus, /j,: 



Ati - M2 = -1.448(si - S2) 



(32) 



Applying Eq. (^ to Coma, ABCG 85, and the group 
of NGC 4839, we obtain: 



A^ABCGSS — MComa 
A'NGC4839 — MC oma 



1.98 ±0.33 
-0.51 ±0.53 



(33) 

The mean entropy are taken from Table ^ while errors 
are deduced from the standard error propagation. 
In term of relative distances we have: 



on 
on 



2.5 
0.8 



+0.4 
0.3 
+0.2 
0.5 



(34) 



The error bars seems especially large in the case of 
NGC 4839/Coma distance. Notice however that the num- 
ber of galaxies used in this case is very low (9 galaxies in 
the group of NGC 4839); the error bars would be ~ 0.2 if 
we had a sample of about a hundred galaxies in the group. 

For comparison, the relative distance between ABCG 85 
and Coma derived from their mean redshifts is ~ 2.4, cf. 
Biviano et al. (1995b), Durret et al. (1998). If we suppose 
that the distance of ABCG 85 is given by (|3^) so that its 
redshift is a composition of Hubble velocity plus its peculiar 
velocity, and a contrario that Coma is at rest in its Hubble 
frame, then ABCG 85 would have a line of sight peculiar 
velocity of ~ -600 ± 240km s"\ 

It seems that this method is not adequate to deter- 
mine individual distances inside a cluster. In Fig. ^we show 
for each cluster the distribution of the distance modulus, 
fi, computed for each galaxy using the Eq. (Bl|) relative to 
the mean specific entropy (the values in Table^. The r.m.s. 
scatter is about 1.1 mag for the rich clusters. Nonetheless, we 
note that using only the a — u correlation, following Young & 
Currie (1995) or Binggeli & Jerjen (1998), and fitting it with 
an arbitrary polynomial of either 3 or 5 free parameters, we 
obtain scatters of the order of 0.83 mag. 

For the Coma cluster we also show the dispersion in rel- 
ative distances (Fig. [ic| ), assuming that the centre of Coma 
is at 137 h~Q Mpc (CoUess & Dunn 1996; Durret et al. 1998). 

The scatter around the Entropic Plane is not produced 
by the distances of individual galaxies with respect to the 
centre of mass, but it is due to a number of effects that will 
be discussed in the next section. This scatter produces an 
artificially elongated cluster; in the present case of Coma, 
we obtain a depth of ~ 1000 Mpc or, neglecting the 5 most 
'distant' galaxies, a depth of ~ 350 Mpc. 



7 DISCUSSION 

From previous work (Young & Currie 1994, 1995; Binggeli 
& Jerjen 1998) and from our own data of elliptical galax- 
ies in two rich clusters and a group, it is clear that there 
are strong correlations among the primary free parameters 




Figure 9. Distribution of the distance modulus, /i, for each clus- 
ter computed using Eq. (pl|). The superposed curves are gaussian 
fits (using a standard least square method) and in each panel is 
quoted the standard error in magnitude. Notice that scales are 
different in each panel. 



of the Sersic profile. Before exploring the possible applica- 
tion of this relationship, we tried to understand the physics 
underlying it. 

Since isolated elliptical galaxies seem to be systems in 
quasi-equilibrium, i.e., only subject to secular evolution, we 
assume that their dynamical state is characterised by a sta- 
tionary entropy. We can go further and postulate that the 
entropy per mass unit (specific entropy, s) would be the 
same for all galaxies. In order to confront this hypothesis 
against observational data, we needed to compute this spe- 
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Figure 10. Individual distances of the galaxies in Coma com- 
puted using Eq. ) and assuming the centre of the cluster is at 
137 h^Q Mpc. Notice how this artificially stretches the cluster in 
the line-of-sight direction. 



cific entropy. From Sersic's light (and mass) distribution for 
elliptical galaxies, and assuming the equations of state of 
an ideal g well as the isotropy of the velocity tensor, 

we have derived a theoretical relationship among the pa- 
rameters describing Sersic's profile. This relation defines an 
'Entropic Plane'. 

We have found that these primary parameters for a 
given cluster or group are well approximated by the pre- 
dicted constant specific entropy relation: elliptical galaxies 
seem to share the same s. 

We would like to emphasise that the two principal as- 
sumptions that we have used, i.e., the entropy and the ideal 
gas equations of state, are not straightforward concepts in 
the context of stellar self-gravitating dynamics. Interest- 
ingly, these supposedly inefficient tools allow us to predict 
relations which match quite well the observational data, a 
result which raises interesting questions on the description 
of self-gravitating systems. 

Nevertheless, we note that the observations are quite 
dispersed along the Entropy Line given by Eq. and that 
there appears to be a tilt between de data and the predicted 
Entropy Line (Fig. A number of effects could account for 
this result. 

• Concerning the model: 

- we have used a constant M/L ratio. This is the sim- 
plest way to account for the relative contribution of dif- 
ferent hght components in a galaxy, and for ellipticals 
as a whole, given our knowledge of the M/L trend in 
these galaxies (j0rgensen 1997). Notice that the specific 
entropy varies as ln(M/L). Such a weak dependence on 
M/L (which we suppose constant from one galaxy to 
another), may certainly explain why the predicted cor- 
relations work so well even if, indeed, there is a variation 
of M/L as a function of v (and consequently of the total 
meiss of the galaxy); 

- for the hydrostatic equilibrium we have considered 
the velocity dispersion tensor to be isotropic. This is 
probably correct for galactic centres, but it is probably 
not for the outskirts. At present, only the introduction of 
a phenomenological parameter describing the variation 



of the anisotropy with galactocentric distance would ac- 
count for such a dispersion; 

- the equations of state that we have used are for 
ideal, non-interacting, hot gas. The inclusion of self- 
gravity explicitly (it is implicitly included by way of 
the hydrostatic equilibrium and Poisson equations) will 
probably tilt the Entropic Plane since the total mass of 
the galaxy correlates with the Sersic parameter, in par- 
ticular i/; 

- galaxies are likely to undergo merging and/or colli- 
sions, especially in clusters; therefore during these events 
they will necessarily escape from the 'Entropic Plane'. 
However, the resulting galaxies are expected to return 
back to the 'Entropic Plane' after a given relaxation 
time. 

• Concerning the data analysis: 

- In Fig. (|l|) we have given the distribution of u; we 
note that this distribution is different from that of Virgo 
galaxies given by Graham & CoUess (1997). On the other 
hand, our distribution of v is similar to the one obtained 
by Binggeli & Jerjen (1998). This seems to indicate that, 
in order to obtain comparable results, it is very impor- 
tant to use the same procedure to analyse the data. 

As noted in the previous section, galaxies in a cluster 
are not exactly at the same distance from us, and the scat- 
ter around the 'Entropy line' is partly due to the various 
(physical and observational) effects noted above. The scat- 
ter is thus larger than that due only to distance variations. 
We have used Formula (^l|) for each galaxy in the Coma 
cluster compared to the mean s (the zero-point for the clus- 
ter). The Coma cluster then appears extremely elongated 
along the line-of-sight; Binggeli & Jerjen (1998) conclude in 
the same vein when discussing the application of the profile 
based distance indicator to individual galaxies in the Virgo 
cluster made by Young & Currie (1995) [but see also Young 
& Currie (1998), who refute the analysis of Binggeli & Jer- 
jen]. 

The 'Entropic Plane' appears to be essentially useful 
to determine the distance to clusters, but not to individual 
galaxies inside a cluster, at least with the present scatter; 
notice that it is also the case for the TuUy-Fisher relation, 
which is used to compare relative distances between clusters 
(Scodeggio et al. 1997). 



8 CONCLUSION 

The distance indicator we propose is not new, in the sense 
of using the distance independent shape parameter (or a 
function depending only of u) of the Sersic law as a yardstick 
(cf. Young & Currie 1994, 1995, who proposed a distance in- 
dicator based on the shape profile of elliptical galaxies and 
applied it to the Fornax and Virgo clusters). Our distance 
indicator, however, was derived from thermodynamics equi- 
librium, based on our strong hypothesis: i.e. the specific en- 
tropy of relaxed galaxies is universal. As noted in the in- 
troduction, this should be a natural consequence of some 
violent relaxation process; the universality, however, seems 
not to be guaranteed since the history or evolutionary stage 
of a given cluster may have some infiuence in the formation 
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process of galaxies and consequently may produce some vari- 
ations in the resulting specific entropy. 

Therefore, seeking a physical explanation for the profile- 
shape distance indicators allowed us to derive a variant dis- 
tance indicator that has only one free parameter, namely 
the specific entropy of galaxies. That should be compared 
with the empirical profile-shape distance indicators [Young 
& Currie 1994, 1995; Binggeh & Jerjen (1998)] that need 3 
to 5 free parameters. Nevertheless, since the exact shape of 
the Entropic Plane depends on the way we model the ellipti- 
cal galaxies, the empirical profile-shape distance indicators 
have a smaller intrinsic scatter. 

Concerning the unique specific entropy of elliptical 
galaxies being independent of the cluster in which they are 
located, if the existence of an universal profile (Navarro et al. 
1995) is confirmed, the peculiar properties of clusters should 
not be very important. In this case, the relation expressed 
by Eq. ( |29[ ) may provide some clue to the relaxation process 
in galaxies in addition to an attractive way to measure rel- 
ative distances between clusters, using only photometrical 
data - in contrast to the Fundamental Plane which requires 
spectroscopic data as well. 
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APPENDIX 1 

In Section we derive the specific entropy from tlie fun- 
damental relation of thermodynamics. Here, show that an 
equivalent form is obtained using the Boltzmann-Gibbs en- 
tropy: 



5^ 



Fln(F)d^2:d^«, 



(35) 



where Q is the volume in phase space and F is the distribu- 
tion function. From the definition of a mean value, Eq. (jssj 
can be re-written as: 



(36) 



The distribution function can also be interpreted as the 
meiss density in phase space, dM/{d^xd^v). Thus it is legit- 
imate to express F as: 

F oc p/a^ , (37) 

where a is the velocity dispersion and p the density. Now, 
we need an equation of state. Adopting the following form: 

P = Kp''a^ , 

(a polytrope of index n = 1/(7 



(38) 



1 



[Binney & Tremaine 
and teq) we 



(39) 



nney 

1987, p. 224]) and substituting it in Eq. ^) 
finally get: 

s = (ln(p-(i+3^/2)p3/2)) . 

For 7 = 1 (ideal gas) we recover Eq. 

APPENDIX 2 



The numerical computation of Eq. (|6|) for an arbitrary choice 
of I/, a and Eq is cumbersome and slow. Here, we will derive a 
simple analytical approximation for the specific entropy for 
a galaxy following the Sersic law. We do this in two steps. 

(I) First we look for a simple functional form that describes 
s as a function of the Sersic parameters. 

In section 2.1 we show that the specific entropy can be 
expressed as a mean value, cf. Eq. m. We rewrite it as: 



5 3 
s ^ ln(pcff ) + - ln(Pcff ) + C , 



(40) 



where we have used the approximation (x) oc Xcff , with Xctt 
meaning the value of x{r) taken at r = RcS, the effective 
radius, and C is a constant of proportionality. Although this 
approximation is justified for homologous, isotropic systems, 
it is in principle not so for the Sersic profile. Our goal here, 
however, is just to obtain a rough analytical form that will 
later be fine tuned. 

Both pcff and Pcff can be computed numerically from 
Eqs. (^, (0) and(|2^). For the effective density, we have to 
a good approximation the following analytical expression: 
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and, for the effective pressure we have: 
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(42) 



Both formulae are analytical fits that give the dependence 
of pcff and Pcff on the parameters of the Sersic profile. 

We also need the relation between the three dimensional 
central mass density, po, and the projected central light den- 
sity, Eq. This can be fairly well approximated as: 



Eo M 0.60 
po ~ p 0.36iy 



(43) 



Substituting Eqs. (^, and (||) into Eq. (^, we 

obtain an analytical expression: 



iln(E„) + |ln(a) + iln(^) +^ln(G) 



+ 



Co + ci \n{u) — — + C2U^ 



(44) 



where the constants c,; are to be determined and fine-tuned. 

(II) The second step is to use the above formula as a tem- 
plate to obtain an accurate approximation independent of 

Pcff. 

This is done by fitting the analytical approximation, 
Eq. (^), to the formula of the specific entropy, Eq. (^ 
(which can be computed numerically). We then get the fol- 
lowing values for the constants c^: 

Co = -0.299, ci = -0.058, ca = 3.656, cg = -1.327. 

The fit is done using a standard least-square method, com- 
puting Eq. (^) in a 'cube' with the axes u, a and Eo . Taking 
the A'l/L ratio as independent of the galactic radius as well 
as constant for all galaxies and scaling G appropriately, we 
obtain the entropy relation expressed in Eq. (Eq), using the 
above numerical values for the Ci's (the numerical value of 
Co has the term (3/2) InG folded in). 

The analytical approximation is accurate to better than 
2% in the range 0.15 < v < 1.85 and the ratios amax/cimin = 
1.2 X 10"* and Eomax/Eomin = 2.3 X 10^ (for a and Eo only 
the ratio is relevant since the values can be scaled with the 
choice of units of G). Considering that the M/L ratio varies 
at most as i\f'''^ (Pahre et al. 1995) and that the range in 
galaxy masses we are interested in is less than 10*, then the 
M/L ratio varies at most by a factor ~ 6. 

Although the functional dependence of the Sersic pa- 
rameters in Eq. (^) was found with the help of the rough 
approximations given by Eqs. (^ and (^), the Eq and a 
dependencies are the same as in the case of a sphere obeying 
the Plummer density profile (where we can calculate exactly 
the specific entropy). 
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